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ABSTRACT 

An analysis is presented for compressible fluid flow across shaft 
face seals and narrow slots. The analysis Includes fluid inertia, viscous 
friction and entrance losses. Subsonic and choked flow conditions can be 
predicted and analyzed. The model is valid for both laminar and turbulent 
flows. Results agree with experiment and with solutions which are more 
limited in applicability. Results show that a parallel film can have a 
positive film stiffness under choked flow conditions. 

INTRODUCTION 

Shaft seal systems in advanced aircraft turbine engines will be 
operated at speeds, temperatures, and pressures higher than shaft seals 
currently used. Conventional face seals presently used in gas turbine 
engines are limited to sliding velocities of about 350 feet per second 
(llO m/sec), pressures of about 125 pounds per square inch (86 N/cm ), and 
gas temperatures of 800 ° F (700 K)(ref. l). Advanced engines, however, 
will require seals to operate to speeds of 500 feet per second (l50 m/sec) 
(ref. 2)5 pressures to 500 pounds per square inch ( 3^0 N/cm ) and tempera- 
tures to 1300'' F (980 K)(ref. 3). For face seals operating at these 
conditions, a positive face separation (no rubbing contact) will be required 
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in order to achieve long life and reliability. A promising approach is 
the use of gas film seals with self-acting lift pads (see fig. 1 and Seal 
Description and Models). 

Since the seals must be pressure balanced, a proper balance of the 
opening and closing forces must be achieved with a leakage gap that has 
tolerable mass leakage. The gap must be small enough so that the leakage 
is minimal but it must be large enough so that power dissipation, due 
to shear in the film, and face deformation are tolerable. Thus the design 
of the sealing gap is vital to seal performance and pressure distribution 
in the gap and mass leakage through the gap must be analyzed. 

In this paper only the sealing dam portion of the seal will be analyzed. 
The classical viscous, isothermal, subsonic, compressible flow analysis 
for this problem is well known (e.g., see Gross, ref. 6). The pressure 
distribution and mass leakage have been calculated for the parallel film 
hydrostatic case. Carothers (ref. 7 ) has conducted compressible flow 
experiments in thin films of air which are flowing radially between 
parallel plates. The pressure distribution was found for both subsonic 
and supersonic axial entrance flows. Carothers' geometry, however, was 
representative of externally pressurized gas bearings rather than seals. 
Grinell (ref. 8) has theoretically and experimentally investigated compres- 
sible flow in thin passages and has shown agreement between theory and 
experiment. Mileller (ref, 9 ) included compressibility for the case where 
there was a restricted by-pass orifice flow into the parallel gap, but his 
analysis did not consider conditions which would yield choked flow in the 
gap. The textbook on mechanical face seals by Mayer (ref. 10) virtually 
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excluded face seal operation under conditions where fluid compressibility 
is important. 

The objective of the report is to present a mathematical analysis of 
compressible fluid flow across shaft face seals and narrow slots. The 
analysis includes fluid inertia^ viscous friction and entrance losses 
and thus describes both subsonic and choked flow conditions. For specified 
pressure ratios , film thicknesses, and friction factor variations with 
Reynolds number, which are used as parametric input, the output variables 
include mass and volume leakage flow rates, force, center of pressure 
and distributions of pressure, density, velocity, temperature and Mach 
number. 


Seal Description and Models 

A conventional face seal which is pressure (force) balanced is shown 
in figure 1. In this seal, the pressure drop occurs across a narrowly 
spaced sealing dam, and the force due to this pressure drop is balanced 
by a hydrostatic closing force and spring force. This configuration, 
however, has an inherent problem. For parallel sealing dam surfaces, 
the force caused by the pressure drop across the sealing dam is independent 
of film thickness for low Mach numbers; hence, there is no way of maintaining 
preselected film thickness which will allow tolerable leakage and still 
have noncontact operation. Since the force is independent of film thick- 
ness, the design also lacks axial film stiffness for sufficient dynamic 
tracking of the stationary nosepiece with the rotating seal seat. In 
order to operate adequately, the seal nosepiece must follow the seal seat 
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surface under all conditions vithout surface contact or excessive increase 
in film thickness, vhich would yield high leakage. Some of these conditions 
are axial runout, misalignment, thermal distortion, coning, and dishing. 

A promising method of maintaining a preselected film thickness and 
achieving axial film stiffness is to add a self-acting gas bearing, such 
as shrouded Rayleigh step lift pads to the conventional rubbing contact 
pressure-balanced face seal (see fig, 2 for a sketch and refs. 4 and 5 
for details of this concept, design and test results). Both the sealing 
dam force due to the pressure drop across the sealing dam and the self- 
acting lift pad (gas bearing) force are balanced by the hydrostatic and 
spring closing forces. The gas bearing has a desirable characteristic 
whereby the force Increases with decreasing film thickness. If the seal 
is perturbed in such a way as to decrease the gap, the additional force 
generated by the gas bearing will open the gap to the original equilibrium 
position. In a similar manner, if the gap becomes larger, the gas bear- 
ing force decreases, and the closing force will cause the seal gap to 
return to the equilibrium position. 

Some Models and Their Limitations 

For subsonic flow, variable cross-sectional area, and Mach number 
<l//p' (all symbols are defined in Appendix A), the analysis of reference 
11 can be used. This differential analysis finds a solution for the case 
when the viscous friction is balanced by the pressure drop and the fluid 
momentum change is negligible. For this condition the entrance pressure 
drop is small, hence ” ^1’ ^2 ~ ^3 flow is 
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nearly isothermal for the film thickness range studied* Neglecting rotation 
and inertia, it can be shown theoretically that the isothermal and adiabatic 
flow solutions on the average are identical (see ref. 12). The analysis 
of reference 11 yields the classical cubic dependence of mass flow on film 
thickness. Also, this analysis enables small tilts of the sealing dam 
surfaces to be studied. These small tilts simulate seal face deformations 
which can occur due to thermal, centrifugal, etc., effects. 

For flow approaching Mach 1, i,e. , choked flow, the inertia force 
neglected in reference 11 becomes important. Also, the flow behaves more 
as an adiabatic flow than an isothermal flow. The equations (partial 
differential equations) are relatively complex because of the nonlinearity 
of the inertia terms; thus an approximate analytical model must be used. 

In the approximate model presented here, the flow is analyzed as a 
quasi-one-dimensional flow utilizing a control volume integral analysis. 

Since in gas film face seal applications there is seal seat rotation 
relative to the seal nosepiece, a question arises as to the validity of 
neglecting the centripetal inertia force which would affect the radial 
flow. This question was examined in reference 13 o It was found that 
the circumferential rotational flow can be uncoupled from the radial 
pressure flow when the ratio of rotational velocity to pressure flow velocity 
is less than l/ J • Rotational effects are important for calculating 
the viscous shear and power dissipation. 

It should be pointed cut that the flow in a gas film face seal is 
qualitatively quite different from the flow in an externally pressurized 
thrust bearing. Both operate with very small film thicknesses, usually 
less than 1 mil. However, gas film face seals are characterized by a 
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small (r^-R^VR^ ratio. Also, the inner cavity is a large supply reservoir. 
Hence, in gas film seals the reservoir pressure does not vary with film 
thickness. On the other hand, externally pressurized gas thrust hearings 
are characterized hy a large (r^-R^)/r^ ratio. (The large surface area 
is necessary for high load capacity. ) In order to maintain positive film 
stiffness, compensating inlet flow restrictors are used. Because of these 
restrictors, the inlet pressure to the hearing varies significantly with 
film thickness. As a result of these differences, the flow regimes for 
the two cases are quite different. See reference 12 for a further dis- 
cussion of this point. 

QUASI -ONE -DIMENSIONAL FLOW ANALYSIS 

As stated in the previous section, when the radial flow is close to 
choked or is choked at the exit, the viscous flow analysis in reference 11 
is no longer valid. Fluid inertia, neglected in that analysis, becomes 
important. The rigorous inclusion of inertial terms necessitates the 
solution of nonlinear partial differential equations and is therefore not 
practical for engineering calculations; hence an approximate analytical 
model will he formulated which will he especially useful for calculating 
flow conditions near and at choking. 

The geometry is shown in figure 3- Two parallel, co-axial, circular 
rings are separated hy a very narrow gap. Rotational effects on the 
radial flow are negligible for most applications (ref. 13). A pressure 
differential exists between the ring's inner and outer radii. The cavities 
on either side (i.d. and o.d.) of the sealing dam are assumed to he 
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constant pressure reservoirs. For sut)sonic flow, reservoir conditions 
and the exit arahient pressure are specified. For choked flow, reservoir 
conditions and an exit Mach number of one are specified. 

The analysis can be separated into two parts, which can then be 
considered separately. One part is an analysis of the seal passage 
itself, in which the flow is assumed to behave as a constant area, 
adiabatic flow with friction, while the other part is an analysis of the 
entrance flow. The analysis for these two regions will be described 
first followed by a discussion of the iterative procedure for the complete 
solution. 


Constant Area Adiabatic Flow With Friction 

It is assumed that the flow in the seal leakage flow region behaves 
as a constant area adiabatic flow with friction. A quasi -one -dimensional 
approximation is made wherein it is assumed that the flow properties can 
be described in terms of their cross-sectional averages. 

The following assumptions have been made in the analysis: 

1. The area expansion due to radius increase is neglected. 

(in most mainshaft face seals the i.do/o.d. is about 0 . 92 L) 

2. The flow is adiabatic. 

3. No shaft work is done on or by the system. 

4. No potential energy is present such as elevation differences, etc. 
The fluid behaves as a perfect gas. 

6. The sealing surfaces are parallel. 

With these assumptions, the flow is commonly known as Fanno line flow 
(ref. 15). This analysis is similar to that used by Grinnell (ref. 8) 


for flow in thin passages. The governing equations when area 
changes are neglected are as follows: 

1^ Conservation of Mass 


[\/^ = = constant 


which reduces to: 


dip , I - n 

p 2 


2, Conservation of Energy 

e + -P + ill = 
? z. 


c CHS tan o 


or 


h 

a 


constant 


where h is the specific enthalpy. This can he written as 

dh= CpdT= 


This equation can he reduced to 
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where M is the Mach number, 
3. Equation of State 

P = f(RT 


( 6 ) 


dP- . dT 

p p r (7) 


4. Conservation of Momentum 

-AdP-ntdAw = Mdll 


( 8 ) 


Nov the following parameters will be introduced: 
hydraulic diameter^ D = 

Darcy friction factor^ f^ = 


It is assumed that the viscous effect can be represented by a mean friction 
factor. The mean friction factor depends on Reynolds number only and 
is a slowly varying function of Reynolds number for many flows, geometries, 
and conditions. 

The set of equations (2, 5^ and 8) can be combined into a single 
equation in terms of the Mach number alone. Details of obtaining this 




equation can be found in Appendix B. The result is 
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Ml M = - (l-M^)clM^ 

This equation can now be integrated from some location X , where 

a 

the Mach number is to the location where the Mach number is 

By the use of partial fractions, the right-hand side of the above 
equation can be integrated to 






dM“= 


( 10 ) 


where 


M 2r i + (11) 

Since the variation of f^ with X is usually unknown, the integral 
of the left-hand side of equation (9) cannot be evaluated exactly. However, 

^ replaced by an appropriate mean value f ^ , the Integration may 



be performed. 
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aU'ix = 


D 




( 12 ) 


The integration of equation (9) therefore yields 


4£:(x,-X0= BfM.)-B(Mt) 

D 


(13) 


In the computer calculation^ it is convenient to define a length 


U = x 


MAX 



such that 

D D 

iHi - BfM„) 

D 


(i4) 


(15) 


(16) 


That is, is the length of the seal (measured from X^^^) for which 

the condition at the exit is sonic (choked). Note that b(M = l) = 0. It 
follovs^ if and are the Mach numbers at the inlet and exit of 

the seal, respectively, that 
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AR= L,-L^ 


(17) 


See figure 4. 

The reservoir pressure and exit ambient static pressure are known. 
It is desirable to find the pressure at any point in the flow leakage 
passage as a function of Mach number. The flow is adiabatic. Thus the 
stagnation temperature is constant everywhere in the seal, i.e., 


(18) 



or (ref. 15, p. 80) 



(19) 


From this relation, the continuity equation, the equation of state, anti 
the definition of the Mach numbers, one obtains 


M,J I + *|lMf 


( 20 ) 
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Entrance Flow 

If it is assumed that the flow in the entrance region is isentropic 
and adiabatic, the analysis is straightforward. The adiabatic energy 
equation is 


h. = h, + 


which can he written as 

x = 1 4 - (y-i)Mf 
t; e 


( 21 ) 


( 22 ) 


By use of the isentropic relation, = constant, and the equation 

of state^ other inlet variables can be found knowing the sealed high 
pressure reservoir values, 



( 23 ) 



{2h) 


ih 


It is ■well known that entrance flows do not behave as isentropic 
flows but that an additional pressure drop is present due to viscous 
friction;^ turning losses, etc* A practical way of accounting for the 
entrance loss is to introduce an empirically determined velocity loss 
coefficient, C (see Appendix C). Hence equation (22) becomes 

Li 



t : 

2C^ 


while equation (23) becomes 


(25) 


R 



zct I 


(26) 


Under the entrance conditions discussed in Appendix C, the entrance 
velocity loss coefficient, C , is related to the head loss coefficient, k, 

Li 

commonly used in hydraulics by the following relation 



(27) 


Iteration Procedure 

The procedure of solving the above equations is an iterative one. 

The solution is found using a computer program which appears in reference 
12. In the iterative procedure, the first step is to solve for the 


minimum film thickness for which the flow is choked. If the actual film 
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thickness is less than this minimum value, the speed of the flow is less 
than sonic at the exit and the exit condition is actual 

film thickness is greater than the minimum for choked flow, the exit 
flow is sonic and the exit condition is = 1, l>ut the exit pressure is 

unknown. The details of the procedure are as follows: 

Choked Film Thickness Evaluation 

When the flow is just choked, it is known that = P^ and M2 = 1 - 
By combining equation (26), the relation for P^/P^, with equation ( 20 ), 
the equation for = (P-L/P3) ^^e seal leakage passage, one obtains 



Since Pq/P^ known, the Mach number M^ can then be determined from 

this equation. Once the Mach number is known, all other flow quantities 
can be determined. The minimum choked film thickness is calculated as 
follows: (1) A film thickness h is assumed (chosen to be 1 mil in the 

computer program). (2) A Reynolds number is calculated from this film 
thickness. (3) This Reynolds number determines the flow regime, laminar 
or turbulent, and therefore the friction factor, (f = const, /Re , where 
the constant and exponent n are program inputs and depend on whether 
the flow is laminar or turbulent, ) (^) A film thickness is then calculated 

by use of equation (16), i.e., 
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4P'AR ^ ZChR 
D h 


(29) 


(For radial flov between co-axial parallel disks and parallel plates ^ 
the hydraulic diameter D = 2h*) If this h does not agree with the 

previous h^ steps 2-4 are repeated until the h^s agree. This defines 

the critical or minimum choked film thickness h*. li should be noted 

that equation ( 28 ) cannot be solved for all values of and jf , 

Consequently 5 h* cannot be found. This condition arises for very low 
subsonic flows; hence the flow is assumed to be subsonic. 

Subsonic Flow Case 

For input film thickness less than the minimum choked film thickness 5 
the flow is subsonic and the exit condition is P^ ” this case it 

is convenient to use the length defined by equation (l4); that is^ 

refer all quantities to a fictitious length of the seal for which the 
condition at the exit is sonic (choked). To start the iteration procedure, 
the entrance Mach number is assumed to be the same as for critical 

flow. All other quantities at the inlet are then calculated. The value 
of the friction parameter B(M^) is found from equation (ll) and the choking 
length from equation (l 6 ), i.e., = B(M^)D/4f^ . The length can be 

calculated from ^2 ^1 " then the friction parameter B(M^) 

from equation (ll). Once B(M 2 ) is known, the Mach number M 2 can be 
determined and, hence, the pressure P 2 . This procedure is then repeated 
with a new Mach number M^ until ~ ^3* 


17 


Choked Flov Case 

For choked flow, the input film thickness is greater than the choking 

film thickness and the exit condition is = 1. As in the subsonic flow 

case, the iteration procedure begins with the entrance Mach number assumed 

to be the choking film thickness entrance Mach number. Again all entrance 

quantities are calculated. The friction factor f’ is calculated from 

^f ^ L 

the Reynolds number and the friction parameter b(m^) = — p calculated 
from equation(ll) . The trial flow length, L, is calculated as 


I _ 6fM.) D 


(30) 


If L is not equal to the true flow length, ^R, the procedure is repeated 
with a new M^ until L does agree with AR. 

Calculation of Variables at Intermediate Points 

Once the entrance conditions are known, the Mach number, temperature, 
velocity, and pressure are determined at each desired location, 
along the seal passage length in the following way: 

1. is found from equation (l4), 

2. ^ 0 ^^) calculated from equation (l6), 

3. M^ is found by an iterative solution of equation (ll), 

4. T^ is found from equation (l9) where T^ = T* and M^ = 1, hence 




( 31 ) 
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5o The velocity, is found from the Mach number definition. 

6, The density is found from the mass conservation equation (l). 



7. The pressure is found from the perfect gas lav, equation (6), 


Computer Program 

The computer program in reference 12 used to carry out this analysis 
is vritten in FORTRAN IV, The program is described in detail including 
all input and program variables and output options. Program listing and 
flow charts are presented and a sample problem with computer printout is 
given. The output can be either the U. S. Customary or International 
system of units , 


RESULTS AND DISCUSSION 

The computer program from reference 12 was used to carry out the 

quasi -one -dimensional flow analysis. By the use of this program, several 

problems have been solved. One such problem is a case used in the design 

study of reference 5. The sealed pressure and temperature were 65 psia 

(45 N/cm^ abs) and 100“ F (311 K), respectively.. The exit ambient pres- 

2 

sure was 15 psia (l0,3 N/cm abs). The sealing dam radial width (seal 
leakage passage lengtl:^ was 50 mils (0.127 cm). Figure 5 shows the calculated 
seal gas leakage as a function of film thickness, which ranges from 0.12 
to 2 rails (0.0003 to 0.0051 cm). 

The analysis agrees with the differential analysis of reference 11, 
which is the classical compressible viscous analysis, for subsonic viscous 
flow when the Mach number is less than l/y/^. As shown in figure 5 , for 
film thicknesses less than 0„3 mil (8 ra), the Mach number is less than 
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1//F = 0.845 and the leakage dependence on film thickness is the classical 
cubic dependence. The isothermal viscous flov model loses its validity 
vhen the Mach number exceeds l/y/^. The present approximate models 
however, is valid for both subsonic and choked flow. Figure 5 shows that 
for gaps larger than 0.3 mil (8>tA.ra) the leakage has a lesser than cubic 
dependence on film thickness. Note that choking occurs at a film thick- 
ness of 0.52 mils (l3/Hm), The limiting case of orifice flow, in which 
the flow rate varies linearly with film thickness, would be achieved 
when the film thickness approaches the order of the sealing dam width of 
50 mils (0.127 cm). 

As indicated in figure 5^ transition from laminar to turbulent flow 
occurs for a film thickness of approximately 1.3 mils (0.0033 cm). A 
Reynolds number, based on hydraulic diameter, equal to 23OO has been 
chosen to be the critical transition Reynolds number. This appears to be 
a universal critical transition Reynolds number for flows in ducts, pipes, 
and bearings. For laminar flow, a mean Fanning friction factor of 
was used. This friction factor is derived from the classical, viscous 
compressible flow solution in reference 12, For turbulent flow, the 
friction factor was chosen to be a constant equal to O.O15O, chosen from 
Egli's data (ref, l4) for annular flow. For the transition flow regime 
the exact nature of the flow, hence friction factor, is complex and not 
fully understood. The friction factor used for this flow regime is 
derived in reference 12 by assuming a smooth transition from the laminar 
friction factor to the turbulent friction factor, Reynolds numbers in the 
range from 2300 to 3OOO are arbitrarily selected as the transition flow 


regime. 
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In order to check the validity of the model for choked flow and for 
a radial geometry representative of face seals, experiments were conducted 
at the Lewis Research Center. 

Leakage flow was studied for two co-axial rings, 5.50 inches (l 4 cm) 
inner diameter, and 6.00 inches (l5.2 cm) outer diameter, separated hy a 
fixed parallel gap of 1.5 rails (O.OO38 cm). The reservoir of 6O.6 psia 

O 

( 4 l .8 n/citi abs) pressure was held fixed and the exit pressure varied. 
Figure 6 shows a comparison of this analysis^ the classical viscous 
subsonic flow differential analysis, and the radial flow experiment. 

Notice that the experiment shows that the flow does become choked. The 
maximum flow rate (choked) is O.OI76 Ibm/sec (O.OO8O kg/sec). The present 
analysis predicts a slightly higher mass flow but agreement is within 19 
percent over most of the range. On the other hand, classical compressible 
viscous flow theory overestimates the flow rate considerably at all pres- 
sure ratios except very near one. The classical theory predicts no 
limiting mass flow. However, if the limit is chosen as the mass flow 
when the Mach number is unity, predicted flows are about 80 percent higher 
than observed experimentally for choked flow. 

The analysis used a friction factor of 24 /^^ 2 h laminar flow 

region. The good agreement in figure 6 indicates that it should be suffi- 
cient for engineering purposes to use this as a predictor. The important 
point is that it was derived analytically from classical compressible 
viscous flow theory. As is the case with duct flow, it appears to be 
nearly invariant with compressibility. Exact agreement could have been 
achieved if a different friction factor had been used. In this case almost 
exact agreement with the experimental curve would be achieved if f were 
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chosen to be 6/Re * ' ; hovever^ it may not be a reliable predictor for 
other experiments and applications. 

There are many possible reasons for the disagreement between the 
analytical and experimental results. It is extremely difficult to maintain 
a uniform film tlickness for the gap sizes of interest as small as 1.5 rails 
(0.0038 cm). Problems in clamping, machining tolerances, and distortions 
due to large pressure differentials result in experimental error. Another 
possible source of disagreement between theory and experiment could be 
due to area expansion, which the theory neglects. In the experiment, the 
inner diameter to outer diameter ratio is 0,92. Thus one could expect 
a lower flow rate than calculated. (Calculated mass flow is based on the 
arithmetic mean diameter.) 

The theory also assumed isentropic entrance flow conditions. For 
large pressure differentials the entrance velocity is large enough to 
cause a substantial entrance pressure drop. The result of accounting for 
this entrance loss condition would also result in a decrease in mass flow. 
As illustrated in figure 6, an entrance velocity loss coefficient equal to 
0.6 gave excellent agreement with experiment. 

The radial flow experiment was repeated for two other reservoir 
pressures. The experimental results are shown in figure 7 and compared 
with the analysis. For a reservoir (sealed) pressure of 40 psia (27-6 

Q 

N/cm abs), the flow was laminar over the entire range studied. Agreement 
with the analysis assuming isentropic entrance conditions was within l8 
percent. With an entrance loss coefficient of 0.^ agreement between 
analysts and experiment was within 5 percent. 
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For the 97*3 psia ( 67 . 1 N/cm abs) reservoir pressure case> the 
transition Reynolds number of 2300 occurred for a pressure ratio of 0.87* 
The turbulent flow regime was set to begin at a Reynolds number of 3 OOO 
which occurred at a pressure ratio of O. 67 . Using the mean friction 
factor of 0.015 from Egli's data yields excellent agreement with the 
analysts. However, for a 0,6 velocity entrance loss coefficient the flow 
was underestimated by 9 percent. Agreement in the laminar flow regime 
was excellent for the analysis with 0,6 entrance loss coefficient. The 
critical transition Reynolds number was calculated to occur at a pressure 
ratio of 0 , 8 . 

2 

Figure 8 shows the 97.3 psla ( 67 . 1 W/cm abs) reservoir pressure case 

experimental results compared with the analysis using a turbulent mean 

0 25 

friction factor of 0, 087 /Re * , The inverse one fourth dependence on 
Reynolds number by the friction factor occurs for many fully developed 
turbulent flows. The analysis with isentropic entrance conditions over- 
estimates the flow by 13 percent; however with an entrance loss coefficient 
of 0.6, the flow was underestimated by only 2.6 percent. In the laminar 
flow regime the curves are the same as in figure 7 ; agreement remains 
excellent. 

Since shaft face seals are pressure balanced, it is necessary to know 
the sealing dam opening force which is found from 


F 


R»-R, 

( P- iX 


(32) 


0 
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For classical viscous, compressible flow, the above equation is easily- 
integrated in closed form and yields 



2W(R^ 

-R.'lR [ 

/P 

1 ^ ] 


3l 

L'-® 

)■] 




Note that the opening force is independent of the fluid properties and 
film thickness. This gives the well known result that parallel surfaces 
yield no film stiffness. However, this is not the case when the classical 
compressible viscous flow theory is no longer valid. Figure 9 shows a 
plot of seal opening force versus film thickness obtained using the computer 
program in reference 12, Notice that for small film thickness (correspond- 
ing to M i//r ), the force is constant, as predicted by classical 
analysis. However as film thickness (and Mach number) increases, the 
force actually increases. This is a condition of negative film stiffness. 

It would be undesirable to operate in this region unless an auxiliary film 
stiffness generating device, such as self-acting lift pads, could make 
the overall seal stiffness positive. This negative stiffness may be one 
of the reasons why conventional gas film face seals "chatter" under certain 
high pressure conditions. 

Also notice in figure 9 that the opening force attains a peak value, 
then sharply decreases with increasing film thickness. This sharp decrease 
indicates a high positive film stiffness which heretofore was not thought 
to be present for parallel surfaces. 

The results using the analysis with an entrance loss coefficient of 
0.6 are also shown in figure 9. The negative stiffness region has been 
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greatly reduced; hovever there is a larger region of positive film stiffness. 
A maximum film stiffness of about 57?800 Ih^/in. (102,000 N/cm) occurs 
at a film thickness of 0,6 mil (0.00152 cm). Notice that there is another 
possible negative stiffness region where the flow is in transition between 
laminar and turbulent flow. The uncertainty of the real flow behavior 
in this region may mean different behavior than that calculated. Positive 
film stiffness also occurs in the turbulent flow regime. Of course, in an 
actual seal design, the advantage gained by this added film stiffness must 
be weighed against the higher leakage flow resulting from choked flow 
operation. 


Limitations of the Analysis 

One of the limitations of this analysis is the neglect of the details 
of the flow in the entrance length. It is common in lubrication theory 
to neglect the entrance region entirely. This is correct for the slow 
viscous flows which characterize lubrication flows. However, when the 
flow becomes choked, the entrance velocities are high and it is possible 
that this neglect is improper. 

The model developed herein may still be applied with good accuracy 
in cases where entrance effects are significant by using an entrance loss 
coefficient (see Eq. (25)-(2?) and Appendix C). Fleming and Sparrow have 
shown in reference l6 the bulk of the entrance losses occur in a small 
region very close to the duct inlet. Laminar incompressible entrance loss 
coefficients have been calculated and measured for a variety of duct shapes 
(for example, refs. l6 and 17). Turbulent loss coefficients are not so 



25 


widely reported; however, they may he calculated for parallel plate 
channels from the Information in reference l8. The values for k are 
generally less than 20 percent of those for laminar flow. Incompressible 
loss coefficients may he adequate for use with compressible flow (see 
Appendix C ) . 

Also the use of the quasi -one-dimensional model can he justified by 
the good agreement with experiments and the short running time on the 
digital computer. 


Seal Design Example 

NASA has designed a seal with self-acting lift pads for potential 
use as a mainshaft seal in advanced gas turbine engines. This seal is 
described in references 4, 5, and 19. The sealing dam portion of this 
mainshaft face seal was designed and studied using the analysis. 

Figure 10 shows the effect of sealing dam radial width on leakage 
flow film thicknesses from 0.1 mil (0.00025 cm) to 1 mil (0.0025 cm). 

The four radial dam widths used in the study were 5 mils (0.0127 cm), 

20 mils (0.051 cm), 50 mils (0.12? cm), and 100 mils (0.254 cm). These 
plots illustrate that, from strictly a leakage point of view, the longest 
leakage path possible is most desirable. However, it is shown in reference 
5 that the leakage path length must be compromised from a force balance 
point of view when surface deformations occur. Also shown in figure 10 
is the leakage flow rate - film thickness relation for a 0.4 mil (O.OOl cm) 
sealing dam radial width. Notice that at film thicknesses greater than 
0,25 mils (0.000635 cm) the leakage coincides with the leakage values 
obtained for a knife edge using the theoretical orifice flow equation. 
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Also notice the linear variation of leakage flow with film thickness 
increase. For this narrow sealing dam width case when flow behaves as 
knife edge flow the velocity entrance loss coefficient is the same as 
the flow discharge coefficient used in orifice flow analysis. 

Figure 11 shows mass leakage flow as a function of pressure ratio 
for a fixed gap of 0,4 mil (O.OOOl cm) and varying sealing dam width. 

Notice the knife edge, as expected, has the largest leakage flow rate. 

These values have been calculated from the theoretical orifice flow 
equation. The computer program result for a sealing dam width of 0.4 mils 
(O.OOOl cm) agreed with the theoretical orifice flow equation. As the 
sealing dam width increases, the critical pressure ratio for choking 
decreases as indicated by the sonic line in figure 11. Again the advantage 
of wider (longer) sealing dams is apparent. 

Recent NASA sponsored tests (ref. 19) of a shaft face seal with 
self-acting lift augmentation have demonstrated the feasibility of operation 
at gas temperatures up to 1200° F (910 k), pressure differentials across 
the seal up to 250 psi (l72 N/cra ), and relative surface speeds up to 
450 ft/sec (130 m/sec). These tests simulated the engine operating 
conditions such as takeoff, climb, cruise, and descent. The seal appears 
to be operating as predicted by the design analysis. Figure 12 shows the 
nosepiece seal assembly after the 120-hour endurance test. Figure 13 
shows a closeup of the carbon nosepiece from the shaft face seal with 
self-acting lift pads after 320 hours of steady state endurance testing. 

The total time of testing on the carbon nosepiece was 338*5 hours. The 
pretesting lapping marks are still observable. A surface profile trace 
indicated that the average wear on the surface was less than 5 micro- 



27 


inches (0.13>W.ni), A few shallow scratches (50 to 100 microinches, 1 to 2 
^m) were noticed. The test was concluded after 500 total hours of carbon 
noseplece operation. The seal faces had encountered over 50 startups and 

”3 

shutdowns. The leakage rates varied from 11 to 32 SCFM (5.19 x 10 to 
15.6 X lO”^ SCMS) with leakage generally averaging about 25 SCFM (l.l8 x 
lO”^ SCMS). These leakage rates were close to those theoretically predicted. 

SUMMAET OF RESULTS 

An analysis has been presented for compressible fluid flow across 
shaft face seals. This quasi -one-dimensional Integral analysis includes 
fluid inertia and entrance losses in addition to viscous friction vhich 
is accounted for "by a mean friction factor. Subsonic and choked flov 
conditions can he predicted and analyzed. The model is valid for both 
laminar and turbulent flovs. The following pertinent results were found: 

1. Results agree with the classical subsonic compressible viscous 
flow theory for Mach numbers less than i/Zr". Excellent agreement with 
experiment is achieved if an entrance velocity loss coefficient of 0.6 
is used in the laminar flow regime. A friction factor of O.OI 5 O in the 
analysis yielded good experimental agreement in the turbulent regime, 

2. Near and at critical flow conditions (Mach number = l)^ the 
analysts shows a negative film stiffness; whereas for choked flow the 
analysis predicts a high positive film stiffness. The analysis with an 
entrance loss coefficient of 0.6 yielded a maximum positive film stiff- 
ness of 57,800 Ib^/in. ( 102,000 N/cra) for a design problem studied. The 
choked flow results contrast to classical compressible viscous flow results 


which show no film stiffness for parallel surfaces. 
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APPENDIX A 
SYMBOLS 


A 

B 

S 

D 

e 

F 

f ^ 

f 

h 

k 

L 

M 

ft 

P 

Q 

r 

AR 

R 

Re 

T 

u 


area 

4f ' L 

friction parameter = — — 

velocity entrance loss coefficient 

specific heat at constant pressure 

hydraulic diameter = 2h 

specific internal energy 

sealing dam force 

Darcy friction factor 

mean Darcy friction factor 

mean Fanning friction factor = f'/4 

film thickness (gap) 

head loss coefficient 

flow length 

Mach number = W/f H T 

mass flow 

pressure 

volume leakage flow rate 
radius 

sealing dam radial width (physical flow length)^ 

gas constant = universal gas constant/molecular weight 

Reynolds number = p u2h/xc 

temperature 

velocity 


flow width 

radial coordinate direction 
absolute (dynamic) viscosity 
specific heat ratio^ 
density 
shear stress 

Subscripts 

N location along flow leakage length 

0 sealed (reservoir) conditions 

1 entrance conditions 

2 exit conditions 

3 ambient sump conditions 
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APPENDIX B 

Friction Factor - Mach Number Relation 

Equation (9) may be derived in the following way. From the equation 
of state (eq. (?)) and the energy equation^ equation 5^ one obtains 



Combining this with the mass conservation equation (eq. (2)) one obtains 

IP „ _ l+(y-i)M^ d,U5 

p a a? (B2) 

Combining this with the momentum equation, equation 8, yields 

Av - [ — n 

V 2 / (B3) 

From the definition of the Mach number, ^ 
equation, one obtains 

( ' - 

Equations (B3) and (b 4) then give 

^ (b5) 


= u' 7m T, and the energy 




(bM 


from vhich equation (9) can be obtained, i.e,, 


4£ldlX = 
D 


dM^ 

ifM'* ( 


APPEiroiX c 


Remark on Loss Coefficient 

There are several effects which cause non-isentropic pressure drops 
in the face seal entrance region. Due to an abrupt geometric change at 
the seal passage entrance there is a non-uniform profile caused by flow 
turning, separated flow, and ''vena contracta" effect which results in an 
entrance pressure drop. The pressure drop in the flow-development length 
is higher than that in the fully developed flow region because of two 
effects (ref. l6). The first effect is higher wall shear caused by higher 
transverse velocity gradients. The second effect is the momentum increase 
as the velocity distribution becomes less uniform. 

The role of the entrance velocity loss coefficient may be better 
understood by considering the incompressible Bernoulli equation relating 
the stagnation reservoir pressure with the static and dynamic pressure 
at the seal entrance. 

^ lOEAL 

The actual velocity at the seal entrance is less than ideal; it can be 
expressed as 


a,-= C^LL, 

IDEAL 

Where C is the entrance velocity loss coefficient (see ref, 20), 
L 

Substitution of equation (C2) into (Cl) and use of the Mach number 


definition yield 
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p . £RMt. 
' 


(C3) 


or 


p = 1? 

1 H- 


(c4) 


A binomial expansion of the denominator of the compressible entrance 
pressure equation (26) yields 



= l + 

2 - 


? 


4- 


4^ C,' 




(C5) 


The first two terms on the right-hand side strongly predominate for 
Mach numbers less than unity. For example, when = 1 and ]p = 1.4, 




r 


+ X(r-i)^“ 

2 



I -|- 6. yO"f~ O.02.-t~ • * ' 


The error is less than 11 percent if only the first two terms are used. 

The accuracy would be greater at lower Mach numbers. The first two terms 
of the expansion are the same as the denominator of the Bernoulli equation 
(c4) for a gas behaving as a quasi -incompressible fluid. Using this 
observation the relationship between the entrance velocity loss coefficient, 
C , and head loss coefficient, k, commonly used in hydraulics can be found. 
The incompressible Bernoulli equation with a head loss coefficient is 


i 
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P_i.£iiL= 

^“ 2 ~2~^ (c6) 

R_(A+i)2|Lf 

=. |3 — ^A-hl ) (C7) 


Comparing equation (C?) vith equation (C3) yields 

Jl=rt “I 

L 


(c8) 


This relation enables one to convert head loss coefficients reported 
in the litera'ture to velocity loss coefficients C , vhich are more conveni- 

li 

ent to use in Mach number relations. 
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Figure 1. - Pressure-balanced face seal with a gas bearing added for axial film stiffness. 



Figure 2. - Pessure-balanced face seal with no axial film stiffness. 
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Figure 5. - Seal gas leakage versus film thickness for parallel sealing dam surfaces. Radial dam width, AR, 50 mils (0. 130 cm); 
sealed air pressure, P^, 65 psia (45 N/cm^ abs); sealed air temperature, Tq, 100° F (310 K); sump pressure, P 3 , 15 psia 
(10 N/cm^ abs). 
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Figure 6. - Comparison of classical viscous flow theory and present analysis with radial 
flow experimental results. 
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Figure 7. - Comparison of present analysis with radial flow experimental results. 
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PRESSURE RATIO, P 3 /P 0 

Figure 8 . - Comparison of present analysis using 0. 087/Re^' turbulent friction factor 
with radial flow experimental results; Pq, 97. 3 psia (67. 1 N/cm^ abs); Tq, 71° F (295 K); 
gap, 1.53° mils (0.0039 cm); 2.75 in. (6.98 cm); R 2 , 3.00 in. (7.62 cm). 
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Figure 9. - Sealing dam opening force plotted as function of film thickness 
for paraliel surfaces comparing classicaf theory with present analysis wim 
both isentropic and 0.6 entrance loss coefficient; Pq, 215 psia (148 N/cm^ 
abs); P 3 , 15psia(10.3N/cm^abs); T^, 800° F (700 K); AR, 50 mils (0. 127 cm). 
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Figure 10. - Effect of seal face radial length on seal leakage; parallel sealing faces; sump pressure P:> 15 psia (10 3 N/cm^ 

abs); sealed pressure, Pq, 215 psia (148 N/cm^ abs); sealed gas temperature, 800° F (700 K). ’ 
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h = 0.4 MIL (0.00102 CM) 

Rj = 3. 265 IN. (8. 293 CM) 
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Figure 11. - Mass leakage as a function of pressure ratio 
for various sealing dam radial lengths. 



Figure 12. - Gas-film seal assembly nosepiece after 120-hour endurance test. 
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Figure 13. - Close-up of carbon nosepiece of the gas-film seal after ?00 hours of 
endurance testing. Total time on seal: 338.5 hours. 
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